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Abstract 

We consider the interlayer pair-tunneling model for layered cuprates, including 
an effective single particle hopping along the c-axis. A phenomenological sup- 
pression of the c-axis hopping matrix element, by the pseudogap in cuprate su- 
perconductors, is incorporated. At optimal doping, quantities characteristic to 
the superconducting state, such as the transition temperature and the supercon- 
ducting gap are calculated. Results from our calculations are consistent with the 
experimental observations with the noteworthy point that, the superconducting 
gap as a function of temperature shows excellent match to the experimental data. 
Predictions within the model, regarding T c variation with interlayer coupling, are 
natural outcomes which could be tested further. 
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1. Introduction 

A salient common characteristic feature of high temperature layered superconduc- 
tors (HTLS) is the existence of well separated Cu02 planes and the presence of strong 
electronic correlations in these planes. This characteristic, together with the observed 
interlayer contact, though very weak, often forms the basis for theoretical modelling of 
HTLS. Many of the to-date models of HTLS are intraplanar, focussing mainly on the 
electornic interaction in a single CuC>2 layer and ignoring, in the first approximation, 
the weak interlayer contact. But unlike these intraplanar models, the interlayer pair 
tunneling (ILPT) model due to Anderson [|IJ is the one where interlayer contact has duly 
been incorporated and interlayer coupling plays an important role in this model to make 
it more compatible with the realistic HTLS materials than other existing models. 

The ILPT model of interacting electrons on coupled superconducting layers could 
naturally yield high transition temperatures (T c ) observed in the layered superconducting 
materials. In the ILPT model, single particle hopping along the c-direction was argued 
to be blocked owing to strong electronic correlations || and the tunneling of only pairs of 
carriers, caused by interlayer coupling, was considered. The idea behind pair tunneling 
is to amplify the pairing mechanism within a given layer, and as a consequence T c gets 
enhanced @. 

It is important to note that, the interlayer hypothesis was based on certain exper- 
imental observations on underdoped cuprates, that is, the c-axis infrared conductivity 
is small, and the c-axis resistivity shows semiconducting behaviour in contrast to the 
a6-plane resistivity which is metallic |4j]. These observations, which are signatures of 
the marginal presence of c-axis transport, led Anderson and co-workers to enforce com- 
plete suppression of interlayer single particle hopping (ISPH) in their ILPT model 0. 
However, later development in sample preparation and the emergence of more accurate 
experimental techniques, made it possible to study the c-axis transport for wide range 
of dopings and it is observed that for sufficiently overdoped cuprates, both c-axis and 
afr-plane resistivities show similar temperature dependence || implying a three dimen- 
sional metallic behaviour of the material. In other words, for overdoped HTLS materials, 
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c-axis transport is metallic-like. This view is further supported by the c-axis optical con- 
ductivity measurements on overdoped cuprates f| |7| , where the presence of Drude peak 
in the spectrum indicates the existence of interlayer single particle hopping. Also, band 
structure calculations yield substantial value for the interlayer single particle hopping 
matrix element for bilayer cuprates ||. 

Thus, guided by the to-date experimental results regarding the c-axis transport in 
variedly doped HTLS, we consider an extended or more generalized version of the ILPT 
model including an effective ISPH between the layers. Consideration of an effective ISPH 
is essential to remain content with the experimentally observed doping dependence of the 
c-axis transport that, the effective ISPH is strongest in the overdoped regime, decreases 
gradually through optimal doping and becomes negligible in the underdoped regime. 
In case of HTLS, the interrelation between the effects of electronic correlations and 
the variation of dopant concentration, also supports this doping dependence of effective 
ISPH. Effects of electronic correlations ||, which inhibits ISPH within the Anderson's 
ILPT model, is strongest in highly underdoped cuprates leading to non-Fermi liquid 
characteristics HlO] , but becomes weaker in the overdoped regime where a Fermi-liquid 



behavior is expected. Thus, it is suggestive theoretically that the restriction on ISPH 
be relaxed in the overdoped regime and ISPH be given due consideration. To remain 
consistent with these observations and arguments, we introduce a probability factor P 
to describe the effective ISPH, such that, P becomes very small or negligible for highly 
underdoped systems and increases with doping leading to significant values of the effective 
ISPH in the overdoped region. 

An important recent development in HTLS is the observation of a normal state pseu- 
dogap. This gap has been inferred from the NMR [II, I2| , optical conductivity |L3], [14 



heat capacity [|TIJ and transport data |RJ in different cuprate materials. In addition, the 
angle resolved photoemission spectroscopy (ARPES) experiments have shown that the 



symmetry of the normal state pseudogap is d-wave like fll7| , similar to that of the su- 
perconducting gap, and its magnitude (E g ) is large for underdoping, but falls off rapidly 
with the increase of doping ]12[ [15|, [l]J. C. C. Homes et ai, by optical conductivity mea- 



surements on Y123 material Jl3|, found that for oxygen- reduced sample (underdoped), 
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where E g is appreciable, the c-axis dc (zero frequency) conductivity decreases with low- 
ering of temperature, whereas at frequencies well above the pseudogap the conductivity 
is temperature independent. Similar evidences are also reported in experiments on other 
cuprates @, [14]| . These findings reveal a clear correspondence between the pseudogap 



magnitude and the suppression of c-axis conductivity. Motivated by these facts we ex- 
press the probability factor P in terms of E g Jl8|, [19| in such a way, that the doping 
dependence of E g makes P to follow the observed doping dependence of c-axis transport, 
i.e. the effective ISPH becomes small for underdoping and significant for overdoping. 
Initially we choose an exponential form for the probability factor P = e~ Eg ^ T , which 
is very small for underdoping (because of large E g ) ensuring a negligibly small value of 
the effective ISPH. With the increase of doping, E g falls off and P increases, yielding 
significant values of effective ISPH for the overdoped systems. Notice that, when ISPH 
probability is finite, a fraction of the available charge carriers will take part in ISPH, 
and only those particles, which are not participating in ISPH, will be available for pair 
tunneling. Since pair tunneling is a two particle process the probabilistic weight for such 
a process is taken as (1 — P) 2 . Here, in the generalised pair tunneling model, we assume 
that the interlayer tunneling of particles could occur via two channels: (i) single partile 
hopping and (ii) pair tunneling, and these two processes are complementary. 

Having introduced the extended ILPT model, we do a mean-field analysis within 
the BCS approximation and obtain self consistent equations for the gap parameter and 
the chemical potential. Focussing mainly on the optimal doping situation (at which T c 
becomes maximum), we calculate T c and the superconducting gap A^ ax (T). Main results 
from our calculations are as follows. Phase diagram of the model (T c versus S) shows that 
the transition temperature T c at optimal doping increases with the increase of E g . T c as 
a function of bare interlayer hopping (t±) shows upward or downward trend depending 
on the values of E g . High values of the ratio of the superconducting gap to T c (~ 6 — 9), 
as observed in HTLS [EtJ, are recovered for realistic range of parameters. Gap variation 
A^ ia:c (T)/A^ ia:r (0) as a function of scaled temperature T/T c shows qualitatively correct 
behaviour as in high-T c cuprates. A rigorous fit to the experimental gap-variation data 



21|] is obtained with a T-linear choice of the probability factor P. 
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The paper is organised as follows. In section-2 we present a detailed account regarding 
the formulation and justification of the extended ILPT model. Section-3 includes a 
brief presentation of the steps involved in the mean-field calculations leading to the gap 
equation. Results from our calculations are discussed in section-4 and section-5 contains 
a summary of results and some comments. 

2. Extended interlayer pair tunneling model: formulation and justification 



The model Hamiltonian for the coupled bilayer complex [19] is given by 



H = E( £ k- - E [V* c^c%c%c§ + h.c 

«,k,er j,k,k' 

" Y,KWc^c%c%c%+h.c}+ E [tfOO^ + H (1) 

i^j,k i^i,k,cr 

which is similar to that of pair tunneling model except the last term that accounts for 
the interlayer single particle hopping. Here, o^J (cj^) is the fermion creation (annihila- 
tion) operator with momentum k and spin sigma, i (= 1, 2) is the layer index, is the 
pairing potential forming Cooper pairs in the afe-plane and fi is the chemical potential. 
The afr-plane band dispersion is taken to be that of Bi2212, obtained from a six param- 
eter tight binding fit [t ,ti,t 2 ,t 3 ,U,t 5 } = [0.131, -0.149, 0.041, -0.013, -0.014, 0.013] eV 
to the ARPES data P2"| . This six parameter band dispersion, used elsewhere |23[, shows 
flat bands in the Brillouin zone and the corresponding density of states (DOS) has a 
power law singularity known as extended van Hove singularity. This is characteristic 
to the high-T c cuprates |24j. Tp ff (k) and tjf(k) represent effective matrix elements for 



pair tunneling and single particle hopping respectively between the layers which in- 
volve the probability factor P. Effective ISPH is taken to be tjf(k) = t b ± (k)P, where 
t b ± (k) = t ± ((cosk x a — cos k y a) /2) 2 is the k-dependent ISPH as predicted by the band 
structure calculations |§ with t± being the bare ISPH matrix element (a is the lattice 
constant). The bare pair tunneling, following the original ILPT model 0, is taken as 
T p (k) = [(ij_(k)) 2 /|ii|] where t\ is the nearest neighbour hopping matrix element of the 
a6-plane band dispersion p2[ , and the effective pair tunneling, as mentioned earlier, is 
represented by T p eff (k) = T p (k) (1 - Pf. 

Regarding the origin of the pseudogap in cuprates no concensus has been reached 



so far. Some authors believe that it is precursor to the supreconducting gap [E9 and it 



represents the pairing energy of preformed pairs without phase coherence. Absence of 
phase coherence among the preformed pairs, is believed to be due to the strong phase 
fluctuations in the quasi two dimensional systems, and superconductivity sets in at the 
phase locking transition temperature of these pairs. However, such a scenario is contra- 



dicted by recent experiments |26], which observe absence of any isotope effect in E g even 
though there remains an isotope effect in T c , suggesting that the interactions responsible 
for superconductivity and the pseudogap are independent and the pseudo gap cannot 
be attributed to short-range superconducting pairing correlations. Similar views are re- 
flected by M. Suzuki et al. by tunneling spectroscopy measurements on Bi2212 material 
p7H - These experiments also found that the normal state pseudogap coexists with the 



superconducting gap below T c . 

A second school of thought pertains to the strong correlation viewpoint for HTLS 
[ ffisfl , where spin and charge degrees of freedom of electrons are decomposed into spinons 
and holons. For underdoped systems, below a temperature Trvb, the spinons are paired 
as Resonating Valence bond (RVB) |29 in the singlet state, and Trvb is much higher 



than the superconducting transition temperature [g(J. The pseudogap may be inter- 
preted as the spin excitation gap over the RVB singlet state. Within this picture, a 
holon has to combine with a spinon to form a real hole which can hop between the 
layers. Cosequently, for underdoped sytems where spinons are paired below Trvb, the 
ISPH is heavily suppressed. Such a picture, which forms the basis of Anderson's pair 
tunneling model, also motivates us to assume that charge carriers, which are available 
for ISPH, decreases with E g . This phenomenon (for the case of short-range RVB) could 
be represented by an exponential probability factor P = exp(—E g /T). Exponential form 
of P signifies that the particles have to be activated above the gap E g to be available for 
ISPH. This activation is favoured as T goes up or E g gets lowered, whereas an increase 
in E g would degrade the activation and lead to the suppression of ISPH. Notice that the 
form of P, as it stands, could mimic the observed doping dependence of c-axis trans- 
port by virtue of the parameter E g , that is, c-axis transport is small for underdoping 
(E g large), gradually increases through optimal doping (E g gets lowered) and becomes 
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significant for overdoping [E g very small). 

The phenomenological behaviour of P, and consequently that of the effective ISPH, 
is also supported by certain existing ideas and results in HTLS. Kumar and Jayannavar 



31J| , while trying to explain the power law dependence of the semiconducting-like c-axis 
resistivity (p c ) of HTLS, showed that the low temperature up-turn of p c could come 
from the T-dependence of the interlayer coupling t± which gets renormalized as a power 



of temperature t e f = t±T a (Eq.(7) of Ref.[[3T]) where a is an exponent of the order of 



unity. This renormalization of t±, which is an adiabatic modification, was originally 
derived by A. J. Leggett et al. considering a coupling of the slow interplanar electron 
tunneling to some bosonic degrees of freedom [j3^] . It has also been established that the 
strong electronic correlations, present in the CuC>2 planes of HTLS in the form of strong 
intraplane electron-electron scattering, could lead to the complete blocking of single 
particle tunneling between Cu02 layers at zero temperature (T = 0), but the tunneling 
of electron pairs remain uninhibited M. These views, if translated in the language of 
our model, would mean that at T = 0, P = and t e f = but ^ 0, whereas at any 
finite temperature (T ^ 0) P ^ and t e f is nonzero. Precisely, the effective ISPH and 
the effective pair tunneling matrix elements within our model show similar temperature 
dependence. 

Thus, our proposed model for HTLS is consistent with the experimental observations 
and theoretical ideas/results regarding the c-axis electrical transport. Effects of intra- 
planar correlations, which is strongest in the too-underdoped region (E g — > oo), can be 
thought of as coming to picture in our model as blocking the ISPH at any temperature 
[[I in the limit of very small doping. As the doping is increased, intraplanar correlations 
become weaker (E g decreases) and ISPH becomes more and more favourable at nonzero 
temperatures, however at T = 0, ISPH remains completely blocked at any doping. Thus, 
in the limit E g — > oo or at T = 0, the probability factor P vanishes leading to t e f = 
and our model becomes equivalent to the original ILPT model, and in this sense the 
proposed model is more general. 

It is important to note, that the presence of ISPH gives rise to bilayer splitting, which 



7 



however has not been observed in Bi2212 system p3| , particularly at low temperatures. 
In this light, we need to check the validity of our proposed model. A study of the 



electronic density of states (DOS) of our model |18| reveals that for optimally doped 
or underdoped systems, when E g > T c , bilayer splitting is absent at low temperatures 
T < 20 K and within an energy resolution of 1 meV. This is because of strong suppression 
of ISPH by Eg. At higher temperatures, broadening of the quasiparticle states would be 
important to obscure such splitting in real systems. A detailed calculation of the ARPES 
intensity curves within our proposed model clearly establishes that bilayer splitting would 
not be observable at low temperatures ~ 13 K as in actual experiment of Ref . |33| . Even at 



a higher temperature ^ 40 K and with a resolution much better than that in experiment, 



splitting remains absent W8 . These outcomes further embolden the consistency of our 



proposed model with experimental results on HTLS. 

3. Superconducting gap equation within the BCS approximation 

Mean field decoupling of the four fermion terms in the Hamiltonian of Eq.(l) yields 

i,k,<x i,k 

+ E [tfMc^c£+h.c] (2) 
with the gap parameter being defined as 

A k = A, k = E V k>k , (c«, iC « ) + if (k) (c%c$) (3) 

k' 

Layer index (i, j) are equivalent, since by symmetry, in-plane pairing is identical in both 
the layers. Presence of interlayer single particle hopping produces two quasiparticle 



bands E k = ^/{e k - fi - tf(k)} 2 + A 2 k and E+ = - // + tf(k)} 2 + A k . 

Self consistent equations, for chemical potential and the superconducting gap are 
obtained as 

i-5 = i-^E(^-^-^)x(^k)-^E(^-^ + ^ ff )x(^k) (4) 

k k 
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and 

E A k'Vi, k , x(JE£))/2 

A k = ^ 1^ (0) 

1-Tf(k)tel+X(^ 



where x(E k ) = IeF^ 11 * 1 ( — 2/' P = ^ n a sca ^ e of kB=l), 5 = 1 — n with n being 
the number of electrons per site and iV is the total number of lattice sites. The pairing 
interaction Vk,k' is separable as V^k' = Vr]^ tj^i, which also makes the k-dependence of A k 
to be shoved in the symmetry factor rj^. Finally, the expression for the superconducting 
gap (from Eq.(5)) becomes 

1 i V l (xiE^) + X (E+)) / 2 

W ~ N k l-Tf(k)( X (E k ) + X (E k + ))/2 

In our calculations, the d x 2_ y 2 symmetry of the pairing state is considered because of 



growing evidence of the same in high-T c cuprates |34], |35[ and this implies the sym- 
metry factor to be 7]^ = (cos k x a — cosk y a)/2. The parameter V, in our case, is 
the nearest neighbour attractive interaction. In experiments, momentum dependence 
of the normal state pseudogap is found to be of d x 2_ y 2 symmetry ]17|], and we take 
E g (k) = E g \ cos k x a — cos k y a\. This introduces a k dependence in the probability factor 
P(k) = exp(-^(k)/T). 

4. Results and discussions 

We numerically solve the self consistent equations (4) and (6) for the superconducting 
gap parameter and the chemical potential. Basically, we focus on the properties charac- 
terizing the superconducting state, since in this state the quasiparticle picture is valid 
at any doping. We calculate the superconducting gap A^ ax (T) and the transition tem- 
perature T c (at which the gap magnitude drops to zero) and also study their variations 
with model parameters. 

4-1. Transition temperature: 

In Fig.l, we plot transition temperature T c as a function of 5, for various values of 
the normal state pseudogap E„ as shown explicitly. For E g = 0, where ISPH takes 
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its full value (P = 1) and the pair tunneling is absent, T c has two peaks. This is a 
consequence of the splitting of DOS due to the presence of unsuppressed ISPH between 



the layers |L8|, [L9|| . As E g increases, the ISPH gets suppressed and the two peak form 
gradually merge into one peak, yielding the maximum transition temperature (T c m ) at 
optimal doing. For realistic values of E g , the one peak form persists which is the case 
in experiments. We choose the bare value of t± = 40meV, in accordance with the 
band structure calculations || suggesting ~ 0.2 — 0.3. Interaction parameter 

is fixed at V = 70meV which makes T c m ~ 100 K for E g ~ T c m . This value of V is 
kept fixed throughout in this communication. Clearly, for fixed values of V and t±, T c m 
increases with the increase of E g and highest value of T c m is obtained in the Anderson 
limit {E g = oo) where ISPH is completely suppressed giving vent to only pair tunneling 
between the layers. It should be noted that E g is kept fixed for every curve in Fig.l, 
whereas, in reality it should change with doping. Here, we are interested in quantities 
at optimal doping (as in Fig. 2 and Fig.3) and hence consider different fixed E g values 
which might represent different cuprate materials. Increase of T c m with E g is consistent 
with the scaling behaviour of the pseudogap ||12| , according to which, at a fixed doping 



the ratio E g /T™ should ideally fall on the same point for different cuprate materials. 

Variation of T c m as a function of t± is shown in Fig. 2. Numerical values of E g for 
different plots are (A, B, C,D,E) = (2, 4, 6, 8, oo) in meV. For small E g , T™ decreases 
with increasing t±. Increase of E g makes interlayer pair tunneling more probable. Con- 
sequently, T c m rises with t± and in the Anderson limit T c m grows rapidly with It 
may be noted that regarding the role of the interlayer coupling on T c , there is mixed 
experimental evidence (see Ref.|J, page 141). Our results show that T c m may increase 
or decrease or may even remain unaffected with the increase of interlayer coupling t±, 
depending upon the values of E g . This is a prediction of our calculations and a system- 
atic study of the variation of T c with t± for differently doped or different class of cuprate 
materials (with changing E g ) is called for. A possible mechanism to vary t± could be 
the application of pressure along c-axis. 
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4-2. Superconducting gap parameter: 

Two important issues relating the superconducting gap, which figures frequently in 
the HTLS literature, are the ratio of the superconducting gap to T c and the temperature 
variation of the gap. Superconducting gap-ratio (gap- width) A^ ax (0)/T c , as a function 
of t± for different E g , is presented in Fig. 3. Alphabetic labels correspond to different E g 
as written in the figure (similar to the values of E g as in Fig.2). In the Anderson limit, 
gap-ratio could reach the value ~ 4.5 for t± within realistic range (~ 30 — 45meV). 
But, for finite E g , the gap-ratio within this range of t± increases as E g is lowered. For 
example, with E g = 4meV one could get the gap-ratio as high a value as ~ 6 for 
t± ~ 45meV. To understand this, note that for any finite E g the probability factor 
P — > as T — > 0. This situation corresponds to the Anderson limit where ISPH is 
completely suppressed. Hence at T = 0, A^ ax (0) is quantitatively same for any E g 
whatsoever for a fixed t±. But for T ^ 0, as E g is lowered, the suppression of ISPH 
gets weaker and correspondingly T c falls off from its value in the Anderson limit which is 
evident from the plot in Fig.2. Consequently, the ratio 2A^ ax (0)/T c is small for E g = oo 
and increases as E g is decreased. 

Finally, We study the temperature variation of the superconducting gap. In Fig.4 
we plot A^ ax (T) / A^ ax (0) as a function of reduced temperature T/T c (dashed lines), for 
different E g values shown. Here t± = AOmeV and the doping is at the optimum level. 
The solid line represents results for BCS superconductors and solid square symbols are 
experimental data |2]J for Bi-cuprate. Clearly, the locus of experimental data is well 
below the BCS-curve. Results from our calculations are all below the BCS-curve, but 
high values of E g makes them closer to the BCS-curve, with the closest one being that 
in the Anderson limit. For small E g our calculations are close to the experimental data. 
However, our results with P = exp(—E g /T) could not match the experimental data for 
the whole range of T/T c . Maximum mismatch is observed in the low temperature region 
where experimental data registers a near linear fall from its T = saturation. 

Anticipating that the linear fall of A^ ax (T) in the vicinity of T = could come from 
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a linear increment of effective ISPH, we consider another possible choice of P as 

T 

P = 7) 

E g + Te~ E !>/ T K) 

Note that the choice of P is not unique. While making a choice, the conditions to be 
satisfied are, P oc T for E g ^> T and P = 1 for E g = 0, which are ensured in Eq.(7). 
Results, with linear-T dependent P and t± = 40meV, E g = 8meV |36 |, are presented 



in the inset of Fig.4 as a dashed line. The matching with experimental data (solid 
squares) is excellent. To our knowledge, this is the first ever calculations that correctly 
reproduces experimental data, affirming the merits of the model under consideration. 
Use of P from Eq.(7) does not qualitatively change other results discussed earlier, rather 
yields higher values of the gap-ratio than those obtained with exponential P for the same 
set of parameters. As for example, t± = 40meV and E g = AmeV yields a gap-ratio 7.1 
for the present case. 

The T-linear form of the probability factor could have phenomenological support 
from the work by Leggett et al, where the renormalization of t± involves a power law 
dependence of temperature as T a with a being of the order of unity. Arguments could 
also be given within the spin-charge separation picture of Anderson and coworkers, where 
spinons and holons are the quasiparticles in a layer, and a holon has to combine with a 
spinon to form a real hole that can hop to another layer. Thus, within this picture, c-axis 
transport is proportional to the spinon density, and since the density of spin excitations 
grows as T within the (long-range) RVB model |29], the c-axis transport is expected to 



increase linearly with T. These views, at least on the phenomenological level, justifies 
the form of P as in Eq.(7). 

So far, we have presented results within the extended ILPT model incorporating the 
pseudogap as a phenomenological parameter which affects only the out-of-plane charge 
transport (suppresses the ISPH), and our focus remained on the optimally doped situa- 
tion (fixed E g ). But in actual case, the pseudogap should also affect the in-plane charge 
dynamics and should vary with doping as observed in experiments jl2[ [15], pjj. Regarding 
the in-plane effects of E g , based on the analysis of experimental pseudogap data for sev- 
eral cuprates, it has been suggested that, the pseudogap suppresses the in-plane spectral 
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weight or single particle DOS [O, and to this effect, the total measurable spectroscopic 
gap on the Fermi surface (FS) would be an additive combination of the superconducting 
gap and the pseudogap. Following this idea, we modify the quasiparticle band enegies as 



= y {ek — jU ± £jf (k)} 2 + A 2 . + E g (k) 2 and do preliminary calculations incorporating 
the variation of E g with doping as in Ref ||12|| . Our results show that, T c in this case, falls 
off rapidly towards underdoping because of the rapid suppression of spectral weight by 
high values of E g . We also find that, the total spectroscopic gap at zero temperature 
(T = 0) shows non-trivial doping dependence (different from T c ) as noted by tunneling 
experiments on HTLS |37|]. Around the optimal doping, in-plane effects of E g does not 
change the qualitative features of the superconducting gap or the transition tempera- 



ture pg| , and the conclusions drawn in this paper remains unaffected. However, to have 
better understanding about the effects of E g on the afr-plane charge dynamics, and to 
draw definite conclusions regarding the behaviour of T c and gap parameter as a function 
of doping within the DOS-suppression picture, detailed calculations would be necessary. 
Results on these matters will be presented in future communications. 

5. Summary and comments 

In this communication, we consider an extension of the interlayer pair tunneling model 
including an effective interlayer single particle hopping and calculate the superconducting 
gap as well as the mean-field T c . Here, we briefly state the results from our calculations. 
Increase of with the increase of E g is in qualitative agreement with the scaling 
behaviour of the pseudogap. Within the model, high values for the gap-ratio are obtained 
and the experimental gap- variation data are reproduced with a T-linear probability factor 
for the effective ISPH. An outcome of our calculations, that T c may increase or decrease 
as a function of interlayer coupling depending on E g , comes as a natural prediction which 
could be put to test. Thus, several contentious issues of high-T c cuprates are explained 
within the extended interlayer pair tunneling model and new predictions are made. 

Certain comments regarding the phase diagram and the extended ILPT model in 
general, would be relevant at this point. In Fig.l we find that for the no doping (5 = 0) 
situation a finite T c is obtained, which is not the case in real layered cuprate materials. 
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This is due to the non-consideration of the effect of the strong on-site repulsions (strong 
correlations) on the kinetic energy term removing any double occupancy. The incon- 
sistency can be removed by calculating T c within a t — J — V model where hopping is 
restricted in a space with no double occupancy [39fl. But the point is that, in both of 



these analyses (in the limit of weak and strong correlations), properties characterising 
the superconducting state would remain unaltered. Moreover, we would like to note 
that, in our extended ILPT model, strong intraplanar correlations do exist, which plays 
role only to suppress the single particle hopping between CuC>2 layers. This is similar in 
spirit to the views adapted by Anderson and collaborators relating to the original ILPT 
model 0. 

In the context of ILPT model, Chakravarty and Anderson showed that, in the limit of 
very small temperature or (T — > 0), the phenomenon of the blocking of ISPH is exhibited 
by those HTLS materials characterized by non- Fermi liquid behaviour fllOfl . We would 
like to note, that non-Fermi liquid behaviour, as mentioned in the introduction, is the 
characteristic of too-underdoped HTLS (E g very large or E g — > oo) and the original 
ILPT model, which is a special case of our extended model, is valid for very little doped 
HTLS. For overdoping, on the other hand, HTLS show Fermi liquid like characteristics 
and the c-axis charge transport becomes more like that in the afr-plane. This implies 
that, for overdoped HTLS, ISPH gains significance and an extension of ILPT model, as 
considered in this communication, becomes essential. Thus one finds that the extended 
ILPT model is applicable to HTLS for a wide range of doping starting from too little to 
very high, whereas the original ILPT model is applicable only in the underdoped region. 
Furthermore, comparison of our results with those from experiments on HTLS show that, 
the extended ILPT model captures HTLS phenomenology better than the original ILPT 
model. 
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Figure captions: 



Fig.l. Mean-field superconducting transition temperature (T c ) as a function of dopant 
concentration (5) for different values of E g (magnitude of the normal state pseudo- 
gap) as shown in the figure. Other parameters are chosen to be V = 70 meV and 
t± = AOmeV. Anderson limit corresponds to E g = oo. 

Fig.2. Plot of T c m as a function of bare interlayer single particle hopping tj_, for different 
Eg values denoted by alphabetic labels (A, B,C, D, E) = (2,4, 6, 8, oo) in meV. 
Varied behaviour of T c m flow with t± is apparent from the curves. 

Fig.3. Maximum value of the gap-ratio (2A^ ax (0)/T c ) is plotted as a function of t± at 
optimal doping. Alphabatic labels corresponding to different E g values are the 
same as in Fig.2 and are listed in the inset. 

Fig.4. Finite temperature gap scaled to its zero temperature value (A^ ax (T)/A^ ax (0)), 
as a function of reduced temperature (T/T c ). Solid line is the BCS-form, solid 



square symbols are experimental data from Ref.[21] and dashed lines are from our 
calculations for different E g shown. [Inset: Dashed line is from our calculations 
with the renormalizing factor as in Eq.(7). Solid squares are experimental data 
and solid line is the BCS-form as in the main figure]. 
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